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Abstract. We present new links between some remarkable martin- 
gales found in the study of the Binary Search Tree or of the bisection 
problem, looking at them on the probability space of a continuous 
time binary branching process. 



1 Introduction 

This paper is a kind of game with martingales around the Binary 
Search Tree (BST) model (see Mahmoud [2H])- Under the random 
permutation model, the BST process is an increasing (in size) se- 
quence of binary trees (T n ) n >o storing data, in such a way that for 
every integer n, T n has n + 1 leaves; the growing from time n to 
time n + 1 occurs by choosing uniformly a leaf and replacing it by an 
internal node with two leaves. In the BST we are interested in the 
profile, i.e. the number of leaves in each generation. A polynomial 
that encodes the profile (called the level polynomial) allows us to de- 
fine a family of martingales (M^ ST (z),n > 0) where z is a positive 
real parameter (see Jabbour 19 J). It is defined in Section ETTT1 

There are (at least) two ways of connecting the BST model to 
branching random walks, and to take advantage of related proba- 
bilistic methods and results. 

The first one consists in embedding the BST process into a conti- 
nuous time process. It is a good way to create independence between 
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disjoint subtrees. The Yule process (see Athreya and Ney 2 ) is a 
continuous time binary branching process in which an ancestor has an 
exponential 1 distributed lifetime and at his death, gives rise to two 
children with independent exponential 1 lifetimes and so on. Define 
the position of an individual as its generation number, and let Zt be 
the sum of Dirac masses of positions of the population living at time t. 
The process (Zt,t > 0) is a continuous time branching random walk; 
call it the Yule-time process. When keeping track of the genealogical 
structure, call TTf the tree at time t and call (TT^i > 0) a Yule-tree 
process. 

By embedding, the BST is a Yule-tree process stopped at r n , the 
first time when n+1 individuals exist (Pittel j2H!> Biggins 7 a , Devroye 
|13|). A continuous time family of martingales (M(t,z),t > 0) is at- 
tached to this model; it is defined in Section 12.21 In a recent paper 
(Chauvin & al. ^2]) where several models are embedded in the proba- 
bility space of the Yule-tree process, the martingale {M^ ST (z),n > 0) 
appears as a projection of the martingale (M(t, z),t > 0). 

In addition, consider the Yule tree and, on each branch, the suc- 
cessive birthdates of descendants of the ancestor. They are sums of 
independent exponentially distributed random variables with mean 1, 
so that it is natural to exchange time and space and to look at these 
birthdates as successive positions in a random walk. Combined with 
the independence between subtrees, it gives a discrete time branching 
random walk; call it the Yule-generation process. 

The second way consists in "approaching" the BST by the so- 
called bisection model (Devroye ^3]) Drmota JED- It is also known 
as Kolmogorov's rock model. An object (rock) is initially of mass 
one. At time 1 it is broken into two rocks with uniform size. At 
time n, each rock (there are 2 n ) is broken independently from the 
other ones into two rocks with uniform size. The mass of each rock 
results from the product of independent uniform random variables. 
Taking logarithms gives a discrete time branching random walk; call 
it the Bisection process. In |12| . it was observed in the Yule-tree 
environment. 

For these reasons, it is worth considering all these models on the 
same probability space. As explained in Section 2, there are three 
branching random walks: Yule-time, Yule-generation and Bisection, 
each one with its family of additive martingales. 

In Section 3, thanks to a convenient adjustment of parameters 
of these four families (the three previous ones and the BST), we 
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establish strong links between these martingales and their limits as 
n (or t) tends to infinity. In Theorems 13.11 and 13.31 we claim that 
there are actually only two different limits (a.s.) in the domain of 
L 1 convergence. On the boundary of this domain, we identify in 
Theorem l3.4l limits of martingales obtained by taking derivatives with 
respect to the parameter. 

To prove these identifications, we need uniqueness arguments 
which are explained in Section 4. We write the (stochastic) equations 
satisfied by the limits of martingales. The solutions of these equations 
have distributions which are fixed points of so-called smoothing trans- 
formations, as defined in Holley and Liggett d8: or in Durrett and 
Liggett ^7j. Owing to known results on uniqueness of their Laplace 
transforms (Liu [21] , [23 , Kyprianou [2J| ; Biggins and Kyprianou [H] ) 
we get equalities in law (Proposition 4.1). 

Section 5 is devoted to proofs of theorems of Section 3. In par- 
ticular, we show that equalities in law are (a.s.) equalities between 
random variables. 

In Section 6, we explain some relations between the above func- 
tional equations satisfied by Laplace transforms and equations stud- 
ied by Drmota in recent papers on the height of the BST ( |141 115j ). 
This allows to get explicit limiting distributions in the case z = 1/4. 

Let us now fix some notation. In the whole paper we are concerned 
with binary trees whose nodes (also called individuals) are labelled 
by the elements of 

U:={0}U \J{0,l} n , 

n>l 

the set of finite words on the alphabet {0, 1} (with as an empty 
word) . For u and v in U, denote by uv the concatenation of the word 
u with the word v (by convention we set, for any u G U, 0u = u). If 
v ^ 0, we say that uv is a descendant of u and u is an ancestor of uv, 
in particular v is the father of vO and vl. We note u >- v to say that 
v is an ancestor of u. A complete binary tree T is a finite subset of 
U such that 

' G T 

< if uv G T then u G T , 
ul GT^uOGT. 

The elements of T are called nodes , and is called the root ; |it|, 
the number of letters in u, is the depth of u (with |0| = 0). Write 
BinTree for the set of complete binary trees. 
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A tree T G BinTree can be described by giving the set dT of its 
leaves, that is, the nodes that are in T but with no descendants in T. 
The nodes of T\dT are called internal nodes. 

2 The four martingales 

2.1 Branching random walks 

A discrete time supercritical branching random walk (in IR) is 
recursively defined as follows: the initial ancestor is at the origin and 
the positions of his children form a point process Z. The distribution 
of this point process Z is a probability on M, the set of locally finite 
sums of Dirac measures. Each child of the ancestor reproduces in the 
same way and each individual also does: the positions of each sibling 
relative to its parent are an independent copy of Z. Let Z n be the 
point process in IR formed by the n-th generation. The intensity of 
Z is the Radon measure [i defined for every nonnegative bounded 
function / by 



and the intensity of Z n is /x* n . We assume 1 < /*t(IR) < +oo (super- 
criticality) . 

We define for 9 e IR 



Let Moo(#) be the a.s. limit. Under a "/clog A;" type condition, we 
have (Biggins' convergence theorem, for instance in 0IBj) 

- if 8A'(9)-A(6) < 0, the convergence is also in L 1 and EMoo(0) = 

1, 

- if 9 A' (9) - A(9) > 0, then M^tf) = a.s. 

By analogy with the Galton- Watson process, we call "supercriti- 
cal" the values of 9 in the first region, and "critical" (resp. "subcri- 





The (positive) martingale associated with this process is 
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tical") if they correspond to equality (resp. strict inequality) in the 
second region. 

In a continuous time branching random walk (in M), the starting 
point is the same as above. Each individual has an independent 
exponential lifetime (of parameter (3), does not move during its life, 
and at its death is replaced by children according to a copy of a point 
process Z exactly as in the discrete-time scheme. The role of A is 
now played by 

HO) = (3(E ( e dx Z(dx) - 1) . 
Ju 

If Zt denotes the random measure of positions of individuals alive at 
time t, the (positive) martingale associated to this process is 

M(t,0)= [ e Bx - tL{ V Z t (dx) . 

Its behavior as t — > oo is similar to the above, with L instead of A 
(Biggins Uchiyama [SOI)- We denote by M(oo,0) its limit. 

At critical values of 9, the limit martingale vanishes. It is then 
classical, since Neveu |27j . to study the family of derivatives J^M n ((9) 
(or JjjM(i, 0)). These are martingales of expectation 0, which con- 
verge a.s. to a random variable of constant sign, of infinite expecta- 
tion, under appropriate conditions (Kyprianou |21j . Liu |25| . Biggins- 
Kyprianou 8 and Bertoin-Rouault 4 ). The details are given below. 

2.2 Examples 

Bisection martingale 

Passing to logarithms in the bisection problem, we get a discrete 
time branching random walk whose reproduction measure is 

Z BIS = 

fi-\o S u + log(i-lT) > where U ~Z^([0, l]) 1 . We have 

A(6) = logE[U- e + (1 - U)- 6 ] = log . 

Let us make a change of parameter, setting z = ^ip, so that A(0) = 
— log z. The corresponding martingale is 

M* IS (z):= eV-^z", (1) 

\u\=n 

1 W([0, 1]) is the uniform distribution on [0, 1], £ (A) is the exponential distri- 
bution of parameter A, ~ means "is distributed as" and =" means equality in 
distribution. 
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where X u denotes the position of individual u. It is easy to see that 
the range of L 1 convergence is z £ (z~,z+), where z~ < z+ are the 
two (positive) solutions of 

2zlogz -2z + l = 0, 

(i.e. d and c in Drmota's notation, see |15j ) 

d c 
z~ = - = 0.186..., z+ = - = 2.155... 
c 2 c 2 

For z = z~ (resp. z+), applying Theorem 2.5 of Liu [25], we see 
that the derivative (M^ IS )'(z) converges a.s. to a limit denoted by 
M^ IS (z) which is positive (resp. negative) and has infinite expecta- 
tion. 

Yule-time martingale 

The Yule-time process is a continuous time branching random 
walk, its reproduction measure is Z = 25_i og 2, and the parameter f3 
of the exponential lifetimes is equal to 1. We have L(9) = 2 1_e — 1. 
The position of individual u at time t is X u (t) = — |u|log2. In- 
troducing the parameter z = 2~ e we have L(9) = 2z — 1 and the 
corresponding martingale becomes 

M(t,z)= 3M e t(i-2*) 

where Zt denotes the set of individuals alive at time t, of cardinality 
Nt- This can be considered as a generalization of the classical Yule 
martingale (e _ *A^, t > 0) which is known to converge a.s. to a random 
variable £ ~ £(!)• The behavior of this family follows the same rule 
as above. Moreover, it was proved in Bertoin-Rouault [3] (see also 
\12\ ) that for z = the derivative M'(t, z) converges a.s. to a limit 
denoted by M'(oo, z), of constant sign and infinite expectation. 

Yule-generation martingale 

The Yule-generation process is a discrete time branching random 
walk and its reproduction measure is Z = 25 e where e l = 8(1) and 
the factor 2 means that the two brothers appear at the same time. 
Since e l == — log U, the intensity n is the same as in the bisection 
case, and then we have A(9) = log again. With the same change 
of parameter, we have a martingale 

M° EN (z)= e^^A (2) 

\u\=n 
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which has the same form as M EIS (z), and has the same range of L 1 
convergence. However the martingales M^ IS (z) and M^ EN (z) do 
not have the same distribution since the dependence between the po- 
sitions X u , \u\ = n, is different in the two models (although the struc- 
ture of random variables along a given branch is the same). Again, 
from Liu [2^j, we have for z = z^ 1 convergence a.s. of (M„ EN ) f (z) to 
a limit M£* (z), of constant sign and infinite expectation. 



2.3 The BST martingale 

A binary search tree (BST) process (for a detailed description, see 
Mahmoud j2H]) is a sequence (T n ,n > 0) of complete binary trees, 
where T n has n internal nodes, which grows by successive insertions of 
data, under the so-called random permutation model. Let us describe 
the dynamics of the sequence of trees. Tree T\ is reduced to the root 
and has two leaves. Tree T n+ i is obtained from T n by replacing one 
of its n + 1 leaves by an internal node and thus creating two new 
leaves. The insertion is done uniformly on the leaves, which means 
with probability l/(n + 1). 

To study the shape of these trees, it is usual to define the profile 
of tree T n by the collection of 

U k (n):=#{u£dT n ,\u\=k} , k>l, 

counting the number of leaves of T n at each level. The profile is coded 
by the level polynomial J2 k Uk(n)z k , for z G [0, oo) and can be stud- 
ied by martingale methods (Jabbour [TJJl, Chauvin & al. [TT1 IT2] ). 
Because of the dynamics of the tree process, this polynomial, renor- 
malized by its expectation, is a ^v„) -martingale, where Tf n \ is the 
cr-field generated by all the events {u E 7y}j< n)Ug u. More precisely, 
we define the BST martingale 



ir __ i)T C n {z) k C n {z) 



where Cq(z) = 1 and 



*w=n£r =<-»■(?)• »>i- 

fc=o v 7 

In the supercritical range z G (z~,z£), this martingale converges 
in L l to a nondegenerate limit (z) and converges a.s. to 
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elsewhere, in particular for the critical values z~ and . For these 
critical values, the derivative (M^ ST )'(z) is a martingale of expec- 
tation 0, which converges a.s. to a random variable M'J^ (z) of 
constant sign, of infinite expectation. 

3 Connections between these martingales 

We now set all these martingales on the same probability space 
and we define below the continuous time tree valued Yule process. 
Roughly speaking, 

- in the BST we keep track of profile, 

- in the Bisection, we keep track of the balance between right 
subtrees and left subtrees, 

- in the Yule-generation, we keep track of time of appearance of 
the different nodes. 

This set-up provides nice connections which are precised in Theo- 
rems O E3] and El 

Let (vt)t>o be a Poisson point process taking values in U with 
intensity measure z/u, the counting measure on U. Let (TT t ) t >o be a 
BinTree valued process starting from TTq = {0} and jumping only 
when (vt)t>o jumps. Let t be a jump time for v.; TTt is obtained from 
TTi_ in the following way: 

if v t £ <9TlV keep TT t = TV and if v t G dTT t _ take TT t = 
TT t _ U {v t 0,v t l}. 

The counting process (N t )t>o defined by 

N t := #8TT t 

is the classical Yule (or binary fission) process (Athreya-Ney t 2J. In 
the following, we refer to the continuous-time tree process (TT t ) t >o as 
the Yule tree process. 

We note = tq < t\ < T2 < . . . the successive jump times (of 
Tf.), 

T n = inf{t : N t = n + 1} . 

We define recursively the time of appearance (or time of satura- 
tion) of nodes by 



S V = Q) SUl ..U n = inf | s > gtil-lh.-! . Vs = Ul „ Un j 
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Figure 1: Continuous time binary branching process. 

(the definition of vt is given above). Actually S Ul " Un is the sum of 
n i.i.d. £ (1) random variables. This yields TL t = {u : S u < t}. 
The natural filtration is (J-t,t > 0) where Tt is generated by all the 
random variables v s ,s < t. Another useful one is {J- n ,n > 1) where 
JF n is generated by the variables S v for all \v\ < n. Finally we will 
use {F( n ) i n > l)i where Tr n \ is generated by TT T1 , . . . , TT Tn . 

This Yule-time process can also be seen as a fragmentation pro- 
cess. We may encode dyadic open subintervals of [0, 1] with elements 
of U: we set 1$ = (0, 1) and for u = u\Ui . . . E U, 

k k 

4=(j>2^2-* + J>2-'-). 

3=1 3=1 

With this coding, the evolution corresponding to the Yule-time pro- 
cess is a very simple example of fragmentation process. This idea 
goes back to Aldous and Shields (PP Section 7f and 7g). In other 
words, for t > 0, F(t) is a finite family of intervals. At time 0, we 
have .F(O) = (0, 1). Identically independent exponential £ (1) random 
variables are associated with intervals of F(t). Each interval in F(t) 
splits into two parts (with same size) independently of each other 
after an exponential time £(1). 

Hence, one has F(0) = (0,1), F(r x ) = ((0, 1/2), (1/2, 1)) where 
t\ ~ £ (1), etc... One can interpret the two fragments I u q and I u \ 
issued from I u as the two children of I u , one being the left (resp. 
right) fragment I u q (resp. I u i), obtaining thus a binary tree structure. 
An interval with length 2~ k corresponds to a leaf at depth k in the 
corresponding tree structure. 



t 
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Figure 2: Fragmentation and its tree representation. 



3.1 Connection Yule-time -> BST 

In ^2]) it is proved that 

(TT Tn ,n>l) '^(T n ,n>l). 

We can now consider the BST process as the Yule process ob- 
served at the (random) splitting times r n . It turns out that M^ ST (z) 
is the projection of M(r n ,z) on J-^ n y It yields nice limit martingale 
connections 



22-1 



M(oo,z) 
M'(po,z) 
where £ ~ £(!)■ 



L-jM* ST (z) for z€(z-,z+), a.s. (3) 



^M£ ST (z) for z = zf, a.s. 



(4) 



3.2 Connection BST — > Bisection 

(v) 

For w6U, let n\ be the number of individuals alive at time t in 
the subtree beginning at node v. In |12| Section 2.5, it is shown that 
the random variables 



Ijiv) 



n 



(«) 



lim ■ 

£ n (father of v) 



(5) 
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are U([0, 1]) distributed, independent along a branch and that U^ v0 ^ + 
jjivi) _ y f or every v. 

It is then clear that we may construct a Bisection branching ran- 
dom walk with these variables. Let us call M^ IS (z) the associated 
martingale: 

M* IS (z):= £ (n^f-V\ (6) 

\u\=n u>v 

Notice that with the notation of Section 2.2, X u = ^2 V , U — \ogU^ v \ 
Theorem 3.1 For z G (z~,z^), and n > 1, a.s. 

E[M EST (z)\^ IS ]=M EIS (z), (7) 

and 

tBIS(^\ _ ji/rBST 



oo 



(z) . (8) 



Theorem 3.2 For z = z^r , we have 



M'* ST {z) ^M^(z). 

It is of course tempting to conjecture equality of the above random 
variables. 



3.3 Connection Yule-time — > Yule-generation 

Let C n := {u : |u| = n} be the set of the nodes in the n-th 
generation of the Yule tree process. For u G U, by definition of the 
Yule-generation process, position X u can be also seen as the time S u 
of appearance of node u, i.e. the sum of the i.i.d. S (1) lifetimes along 
the branch from the root to u. So 

MZ EN (z):= e {l - 2z)SU z n . (9) 

The following theorem is analogous to Theorem l3.11 it is valid in the 
supercritical case. Theorem 13.41 concerns the critical case. 

Theorem 3.3 For z G {z~ ,z^), a.s. 

E[M(^,z)\r n ]=M^ EN (z), (10) 

and consequently 

Mg EN (z) = M(oo,z). (11) 
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Theorem 3.4 For z = z^r , a.s. 

M'£ EN (z) = M'(oo,z). (12) 

Remark 3.5 

A set L of nodes is usually said to have the line property if no node 
of L is an ancestor of another node of L. In other words, the subtrees 
starting from nodes of L are disjoint trees. Theorems 13.31 and 13.41 
could appear as a consequence, for the particular lines Z t and C n , 
of a more general theorem which would be: additive martingales 
associated to a sequence of "lines" tending to infinity have the same 
limit, independently of the choice of this sequence. This theorem 
holds without any serious difficulty as soon as the notion of "line" 
is precisely defined, which is necessary since several notions exist 2 ; 
"optional lines" in Jagers [20] require measurability of the stopping 
rule with respect to the process until the line. More restrictively, 
"stopping lines" in Chauvin |10| and Kyprianou [221 123j or frosts (in 
the fragmentation frame, cf Bertoin require measurability of the 
stopping rule with respect to the branch from the root to some node 
of the line. The above mentioned general theorem holds for stopping 
lines and not for optional lines in the Jagers' sense. Of course Z% and 
C n are stopping lines, but the stopping time r n (the first time when 
n intervals exist in the fragmentation) defines an optional but not 
stopping line. 

In other words, in view of equalities M EST (z) = M EIS (z) and 
M(oo,z) = Mg EN {z) and in view of connections (J3J) it is clear that 
(for z / 1/2) M(oo,z) is different from M EST (z). This is consistent 
with the fact that r n does not define a stopping line. 

4 Smoothing transformations and limit dis- 
tributions 

Before proving almost sure equalities announced in the theorems 
of Section 3, we first look at equalities in distribution. The random 
limits mentioned in Section 2 satisfy "duplication" relations, which 
come from the binary branching structure of the underlying processes. 
These relations may be viewed as equalities between random vari- 
ables or as functional equations on their Laplace transforms. The 

2 let us also mention the close notion of "cutset" in Peres |28] 
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corresponding distributions are fixed points of so-called smoothing 
transformations (Holley and Liggett ^Hj> Durrett and Liggett [TTj 1 ). 

4.1 Duplication relations 

1) Let us first consider M(oo,z). After conditioning on the first 
splitting time t\ of the Yule-tree process (recall that t\ ~ £(1)), we 
get (HU Section 3.1), for z G (z~,z+) 



M(oo, z) = ze {1 ~ 2z)Tl (M (oo, z) + Mi(oo, z)) a.s., (13) 



where the random variables Mq(oo, z) and Mi(oo, z) are independent, 
distributed as M(oo, z) and independent of t\. Moreover P(M(oo, z) > 
0) = 1. Iterating fTjft we get 



For z = z^r, M(oo, z) = a.s. and the relation satisfied by M'(oo, z) 
is the same as (|T3*|) (mutatis mutandis). Moreover P(M'(oo, z~) > 



2) By definition of the Yule-generation process @, we have, con- 
ditioning upon the first generation, 



Mg EN (z) = (M^(z) + M?S?{z)) a.s., (15) 



which is exactly the same equation as (|13|) . The same result holds 
for derivatives at z = zf. 

3) Let us see now what happens for the BST martingale limit. By 
embedding, it is shown in Section 3.1 that 



Mg ST (z) = z (u^M^z) + (1 - Uf^M^z)) a.s., (16) 



where U ~ W([0, 1]) is nothing but as defined in (jSJ), where 
M^f^(z), M^fZJz) are independent (and independent of U) and 
distributed as (z). For z = zf the relation is the same with 




(14) 



0) =P((M'(oo,z+) < 0) = 1. 



M^ ST (z) instead of Mg ST {z). 
Iterating (|lfjjl. we get 



m° st {z) = z« y: ( n c/M ) 2 



:z-i 




(z) a.s. 



(17) 
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4) By definition of the bisection @ 

m* is (z)= y, {U u(v) ) 2z ~ lzn a - s - 

\u\=n v^u 

We have, conditioning upon the first generation of this process, 

Mg IS {z) = z (u^M^z) + (1 - Uf^M^z)) a.s. (18) 

which is exactly the same equation as (|16|). 

At this stage, we see two packages (M(oo, z), M^ EN (z)) and 
(M^ ST (z), M^ IS (z)), which are consistent with the results of Theo- 
rems an d E3J but which do not yet give equality in law, owing to 
lack of uniqueness. 

4.2 Functional equations and equalities in law 

The relations between random variables in the previous subsec- 
tion imply distributional equations which can be viewed as functional 
equations on their Laplace transforms. 

Set 



3 


(*. 


x) 


: = Eexp 


— xM(oo, z) , 


for z G {z~ 


z + 


3 


(*, 


x) 


:= Eexp 


-x\M'(oo,z)\ , 


for z = z c , 




jGEN 


(* 


x) 


:= Eexp 


-xMg E »{z) , 


for z G (z~ 


Z + 


jGEN 


(z 


x) 


:= Eexp 


-AM'£ EN {z)\ 


, for z = z c 




3 BST 


{z 


x) 


:= Eexp 


-*Mg ST {z), 


for z G (z~ 


Z + 


3 BIS 


(z 


x) 


:= Eexp 


-xM* IS (z), 


for z G (z~ . 


4 



Let us summarize some results on hxed points of smoothing trans- 
formations that are needed for our study. There is a broad literature 
on this topic. One of the more recent contributions is in Biggins and 
Kyprianou [B] • We choose to give these results under the assumptions 
of Liu [2U, US] (see also Kyprianou jUj), which are fulfilled in our 
examples. 

Let us consider a branching random walk as in Section 2.1 with 
Z = Y2i=i $Xi an d P(N = 0) = 0. We defined the martingale 

M n (0) = [ e dx - nA ^Z n (dx) , 
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and the derivative 

M' n (9) = f (x-nA'(e))e ex - nA ^Z n (dx). 

Let us assume EN 1+S < oo and E[Mi(#)] 1+5 < oo, for some 5 > 0. 

In the supercritical range, i.e. if 9A'(6) — A(9) < 0, the Laplace 
transform J(s) = Ee _sM °°^ satisfies 

N 

J(s) =EjJj(se eXl - A(e) ) (19) 

i=l 

(branching property at the first splitting time) and 

Jim 1 ~ J(X) = 1 (20) 

(L 1 convergence of the martingale). 

Moreover for every K > there is only one solution of (|19j) in the 
class of Laplace transforms of nonnegative (non degenerate) random 
variables satisfying 

Jim i-Jfo) = Km (21) 

^■^0+ x 

For critical values, i.e. for such that 9h!{6) — A(9) = 0, the a.s. 
limit M^ )O (0) is positive if 6 < and negative if 9 > 0. Its Laplace 
transform J(s) = Ee~ s ' M °°' 9 " satisfies 

N 

J(a) =EjJj(se te< - A W), (22) 

i=l 

and 

Hm Lzj[M = |g|-i ; (23) 
x^o+ ac| log x\ 

(see Theorem 2.5 a) of Liu |24| with a slight change of notation). 
Moreover for every K > there is only one solution of (|22j) in the 
class of Laplace transforms of nonnegative (non degenerate) random 
variables satisfying 

lim 1 ~ 3 ( X ] = K. (24) 

x^0+ X\ log X\ 

In our setting, this yields the following identities. 
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Proposition 4.1 a) For z G (z c ,z£), we have j(z, •) = j GEN (z, •) 
or equivalently 

M(po,z) l = Mg EN {z) . (25) 

b) For z G (z~,z+) ; we foaue j BST (z,-) = j B1 (z, ■) or equiva- 
lently, 

Mg ST {z) l ^Mg IS {z). (26) 
For critical z we have 

Um w7(f.*) =Um w;"w.*) = 2 

x|0 x|logx| xj.0 x|logx| 2^c — 1 

Um Lzj^f! = lim l-Jf"K-.'> = . (28 , 

xj,o sc|logx| xj.o s | log a; | 1 — 2z c 

Proof: Recall (see end of Section 2.2) that the branching random 
walks BIS and GEN are different but have the same A(0), giving then 
two versions of (|22|). They share the same critical points (recall the 
correspondence = 1 — 2z). 

a) From 1|13|) and (|15jl it is clear that for z E [z~ ,z£], the functions 
j(z, •) and j (z, •) are non-constant solutions of 

J(x) = / J(zxu 22 - 1 ) 2 du, (29) 
J o 

J(0) = 1. (30) 

Moreover, since EM(oo,z) = EM GEN (z) = 1 for z supercritical, it 
turns out that they satisfy (|2U|) and then they are equal. 

b) From (Q2J) and (0 it is now clear that for z G [z c , z c 1, the 
functions j BST (z, •) and j B (z, •) are non-constant solutions of 

J(x) = J 1 3 {xzu 2 *- 1 ) 3 {xz(l - u) 2z - 1 ) du , 

J(0) = 1° (31) 

With the same remark as above, we have uniqueness. 

c) We see that (J23J) yields (|27 |) and ffl ty . ■ 



To end this section let us notice that for the Yule-tree and the 
BST, the limit martingale connection (j3J) gives an important relation 
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between the Laplace transforms: 

J(z, x) = j BST (z, e^dr, . (32) 

5 Proof of theorems 

5.1 Proof of Theorem 13.11 

The relation follows from |JTT|). ©, independence of M BS J(z) 
with respect to (U^ v \v -< u) and the fact that EM^J(z) = 1 for 
every u £ U. 

To prove (JSJ we first pass to the limit in n to get a.s. 

m^)\^ ia \=Mg I8 {z). (33) 

Set for a while, X := Mg ST (z), Y := Mg IS {z) and A := F^IS _ 
Summarizing (|33() and (|26|). we have 

E[X|^] = y and I ! = 7. 
From Exercise 1.11 in [9: this implies X = Y a.s. ■ 



5.2 Proof of Theorem 1331 

It is exactly the same line of argument as in the above subsection, 
using (US]) instead of (JUJ) and (O instead of (fT7|). 

5.3 Proof of Theorem HOI 

Since we work with fixed z, we omit it each time there is no 
possible confusion. 

The idea is to take advantage of the L 1 convergence of a mul- 
tiplicative martingale and then come back taking logarithms. The 
proof can be given both in supercritical and in critical cases and we 
choose to present it for both cases, because it is not more compli- 
cated. Thus it will give an alternative proof of (fTT|) and the proof of 

• First step: multiplicative martingales. Notice that it could be more 
or less directly imported from theorem 3 in Kyprianou |22| but ne- 
vertheless we give details to make the proof self-contained. 
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Multiplicative martingales have appeared many times in the lite- 
rature, for instance in Neveu hi Chauvin jlUj in the branching 
Brownian motion framework, in Biggins and Kyprianou [S] for dis- 
crete branching random walks and in Kyprianou (22j for branching 
Levy processes. They are studied in their own right in relation with 
functional equations or smoothing transformations and also, like here, 
to help find information about additive martingales. Recall that here 

M% EN {z) = Y, ^e (1 - 2z)S \ 
M(t,z) = Y z M e t(1 - 2z) . 

u£Z t 

Let for any real y, 

V n {y)= I]j(^ s "), 

\u\=n 

where j(x) = j GEN (z, x) (it is a solution of equation (|29j) with initial 
condition PU|0 . and let 

r(t)(y)= n jivzMel 1 -**). 

u£Z t 

To prove that V(t)(y) (resp. V n {y)) is a Tt (resp. J rn )-martingale, 
decompose the set Z t (resp. C n ) with respect to a preceding line 
Z s (resp. Cn-i), apply the branching property and get the property 
of constant expectation: decompose E(P(i)(y)) (resp. K(V n (y))) ac- 
cording to the first splitting time t\ and use the fact that j is a 
solution of equation (|29|). 

bmce < j < 1, the martingale (V(t)(y),t > 0) converges when 
t goes to infinity a.s. and in L 1 to a limit V(oo)(y); the martingale 
{j- > n{y)i'n > 0) converges when n goes to infinity a.s. and in L 1 to a 
limit Poo(y). 

Let us now see why these two limits are equal: divide the set of 
individuals alive at time t into those whose generation number is less 
than n and those whose generation number is greater or equal to n. 
This gives a decomposition of V(t)(y); condition with respect to T n 
and apply the branching property; fix n and let t tend to infinity to 
get 

V n (y)=E(V(^)(y)\r i ) a.s., 



Connecting tree models via martingales 19 

and finally, since V(oo)(y) is J 700 := V n .F" measurable 

V 00 (y)=V(oo)(y) a.s. (34) 

• Second step: back to the additive martingale, taking logarithms. 

Let (t n ) n >o be a sequence going to infinity when n goes to infinity. 
We use the behavior of j(z, .) near 0. 

Recall that for z fixed, when i { 0, Mmj(z,x) = 1, so that 
— log j(z,x) ~ 1 — j(z,x). Moreover from (j2U|) and l(T7|) we have 
the sharp estimates 

lim — = 1 for z € (z~, 2d") , 

x—>0 x 

l-j(z,x) ± 
lim — — = Kq tor z = z~r , 

£— »o a; I log a; I 

where Kq = 2/\2z — 1|. The quantity 

m n (z) := m&x{zMe tn{1 - 2z) ; u G Z t J 

satisfies m n (z) < M(t n ,z). If z = z^, taking into account that 
lim n M(t n , zf) = 0, we have \im. n m n {zf) = a.s. Now, for z super- 
critical we check easily that 

m n (z) < m„(z+) logz/logZc+ if 1< z < z+ , 
m n (z) < m n (z c -") log2 / los2c " if z~ < z < 1 , 

and then \\m n m n (z) = a.s. Everything holds in the same way for 
the n-th generation and for M^ EN (z) instead of M(t n , z). 

We deduce that for every e > there is some no such that for 
every n > uq and u £ Z tn (resp. C n ) 

y(l- e ) z M e Mi-2 2 ) < -\ og j(z,yz^e^- 2z ^ tn ) < y(l + e ) z M e Mi-2*) 
in the supercritical case and 

(l-e)y(|n|logz + t n (l - 2z)) Z HM X ~^ 

+ Ko(l-e)y(logy) 2 H e ^(i-2^) 

< -logj(z,i/«M e (l-2*)t») 

< (1 + e)y(\u\ log z + t n (l - 2z))zl"l e *™( 1 - 22 ) 
+ K (l + e)y(logy)z^e t ^ 1 - 2 ^ 
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in the critical case (resp. the same with S u instead of t n ). Adding 
up in u we get respectively 

yM{t n ,z)(l-e) < -\ogV(t n ,y) < yM(t n ,z)(l + e) 

and 

yM'(t n ,z)(l-e) + K (y log y) M(t n ,z) (I + e) 

< -logV(t n ,y) 

< yM'(t n , z){\ + e) + K (y log y)M(t n , z) 

(resp. the same - mutatis mutandis - for the GEN). Taking limits in 
n, this implies a.s. 

- log V(oo, y) = yM(oo, z) , 
-log Poo (y) = yMg EN (z), 

for z supercritical and, for z critical, the following 
-logP(oo,y) = yM'(oo,z) , 

-logPoo(y) = yM'£ EN (z). 

With (|H4|1. we now may conclude the proof. ■ 

5.4 Proof of Theorem EO 

Let us consider only z = z^ to simplify. We know that j B (z, •) 
and j BIS (z, •) satisfy the same equation Q31JI . Since we know by (|2*7|) 
that j BIS (z, •) satisfies 

x^o+ x| log a: | 2z? - 1 

it is enough to prove that j BST (z,-) satisfy also (uniqueness 
mentioned in Section 4.2). By Theorem 13.41 and l)27|) again, j(z,-) 
satisfies also (|35|). Now j B (z,-) is connected to j(z,-) by ((3*2l . 
From some elementary calculations and known properties of Laplace 
transforms we conclude that j BST (z, •) satisfies also ([35 j). ■ 



Connecting tree models via martingales 



21 



6 Links with Drmota's equations 

In this section, we make precise some probabilistic counterparts 
of solutions of equations introduced by M. Drmota in |15L 1161 114) as 
an analytical tool for a sharp study of the height of BST. 



6.1 General case 

We need some changes of parameter, variables and functions. For 
z G [z~ ,z+\ with z 7^ 1/2, set 



a(z) = z 2 



1 

z-1 



For z~ < z < 1/2, (i.e. c' < 2z < 1) the function a increases from 
a c > = e 1//c to +oo, and when 1/2 < z < z+ (i.e. 1 < 2z < c) it 
increases from to a c = e 1//c . We often write a instead of a(z) to 
simplify. 

For 2z l let ip(z,x) = x~ 1 j(z,x 1 ~ 2z ). Equation (|29j) is trans- 
lated into 

/•oo 

ip(z,x) = a~ 2 / (p(z,y/a) 2 dy . (36) 

J X 

Moreover, equation l|3()j) becomes: 

lim x(p(z, x) = 1 , if 2z > 1 , (37) 

x— >oo 

lim xtp(z, x) = 1 , if 2z < 1 . (38) 

x— +0 

and equation (|2()j) becomes: 

lim 1 ~ x ^ x ) = i , if l<2z<c, (39) 

x^oo x 

l im 1 ~M z ,x) = c ' < 2z < 1 . (40) 

Drmota used the solution of the retarded differential equation 

1 



$'(s) = ? $(x/a) 2 , (41) 



a 2 

$(0) = 1, (42) 

and also the solution of the (retarded) convolution equation 

r-y 

y^f(y/a) = / ^(w)^(y - w)dw . (43) 
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If tp satisfies (|1T|| . then <p K denned by ip K (x) = Kip(nx) satisfies the 
same equation. Similarly, if ip satisfies (j43|) . then ip K defined by 
ip K {u) = i{j(u/k) satisfies the same equation. Drmota (Lemmas 18, 
19 and 23 of [IE] and Prop. 5.1 of [13]) proved that for 1 < a < a c , 
there is a unique entire solution (p a of (|4*T|) . and that (p a is a Laplace 
transform: 

/>oo 

<p a (x)= / My)e~ xy dy, (44) 

JO 

where ip a is solution of (|43jl. Moreover, these functions have the 
following behavior: 

Ki, if 1 < a < a c , (45) 
K2 , if a = a c , (46) 

where K\ and £ (0, 00). This implies in particular that </? a satisfies 
(|36f) . But, if <p satisfies (|1T|) . then i/? K defined by (£> K (x) = mp(nx) 
satisfies the same equation, so that the role of the initial condition 
1)42(1 is to fix the constant k. Conversely, we may choose k(z) such 
that V 9 ^) satisfies (|43|) (resp. (|4f)|0 with limit 1 instead of K\ (resp. 
K2). Comparing with (|39[) . we have (at least for 1 < z < z£) that 

k(z) / \ 

For 1 < z < Zc , we know from (|25|) and Theorem 13.31 that the 
Yule-time limit martingale M(oo,z) and the Yule-generation limit 
martingale M^ EN (z) have the same Laplace transform. We see now 
that 

j{z,x) =j GEN (z,x) = X 1^21 K (f a (KX 1^21) , (47) 

or equivalently 

Ee -^M(oo,,) = Ee -^Mg E N {z) = Kx(pa{Kx) (48) 

for some constant k > 0. 

In the critical case, z = z+, we use the derivative martingales. 
From Theorem 13.41 and Section 4.2 we see that (|47|) holds again or, 
in other words, 

Ee _,i- 2 *eV'(oo,* c +)| = Ee -*W|^*"(,+)| = KXipac{Kx) . (49) 



lim 

X — >oc 

lim 



x— >oo x 



1 — X(/? a (x) 
1 — X(/5 a (x) 

log X 



l-2z 
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for some k. 

Now let us consider ijj a as defined in (|l3*|) or (|44|). The relation 
(|32|) together with (jSJ) and (JUJ) gives easily 

*.(,/«) 

The BST limit martingale M^ ST (z) and the Bisection limit martin- 
gale M^ IS (z) (which are equal in the supercritical case) have the 
same Laplace transform and 

In the critical case, z = z£, we have (recall Theorem I3.2j) 

-y^lzL\M'BST (z + )l _£!£-i| M /s/s ( + } | 
Ee r ( 2,+) °° ' c =Ee r (2 Z +) °° = Tp a ^ y j K y 



6.2 Particular cases 

1) For z = 1/2, we have M{t, 1/2) = 1 and j(l/2,x) = e~ x . In 
this case a is not defined. 

2) For z = 1, as previously mentioned, we have M(i, 1) = e~ f Nt 
whose limit is £ ~ £(1), of Laplace transform j(l,x) = ]4^- This 
yields ip(l,x) = jt^., which corresponds of course to ip a (x) with 
a = a(l) = 1. 

3) The function $(x) = e~ x ^ is solution of (JHJ) and (g2I) for 
a = 2, but it does not correspond to the limit of a martingale from 
our families. In fact, there is no real z G (0, oo) such that a(z) = 2. 

4) Drmota ^H] m his Lemma 11 noticed that, for a = 16 

satisfies equation (j41|) and the initial condition (|38|) . 

In (0,oo) the equation q(z) = 16 has two solutions: z = 1/8 and 
z = 1/4. Since 1/8 < z~ , the corresponding martingale M(t, 1/8) 
goes to a.s. and then ip{l/&,x) = 1. Besides, we have z~ < 1/4 < 
1/2, so that j(l/4,x) is not constant and satisfies (J2HJ), (GIB, (HH- 
Again, for every constant k, the function (p K = K(p{n.) satisfies the 
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same system. The constraint (|4U[) leads to k = 4. Taking into account 
the uniqueness, we see that 



Eexp-xM(oo,l/4) = (1 + V2x)e 

which allows identification of the law of M(oo, 1/4). As it is noticed in 
another context in Yor pp. 110-111, this is the Laplace transform 
of the density 

G(t) = —L= e~ l / 2t , t >0, 

of (273/2) -1 ) where we use the notation 7 Q for a variable with distri- 
bution gamma with parameter a, i.e. 

P( 7q G <%) = J-^y^e^dy y > . 
T(a) 

In other words 

M(oo,l/4)^(2 73/2 )- 1 . (50) 

Let us now look for the law of M^f T (l/A). It is easier to work with 
A := (M^ ST (l/4))~ 2 . The limit connection © yields 

(273/2) 2 l = TTfl'A, 

where A is independent of 71 (which is our £). Introducing the Mellin 
transform we have 

E(273/ 2 ) 2S = ^E( 7l ) s E(^). (51) 
Recall that E( 7a ) s = for Rs > -a, so that, for "Sis > -3/4, 

_ 2 2 -r( 2s + 1) 1 
E(A) ~ r(§) 7r*r( a + i)' (52) 

With the duplication formula 

r(2y) = (2^)- 1 / 2 2 2 ^5r(y)r(y + ~) (53) 
taken at y = s + § and y = | , the identity Q52JI becomes 

EUn " r i6 V r(g + !) r(S + |) (54) 

( )_ UJ r(|)r(, + i) r(|) " (54) 
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Now, for a, b > 0, a beta random variable (3 a f, with density 
F(^ b Gdy) = ^±^-y a -\l-y) b ~ 1 , y € (0, 1) , 

r(a)r(&) 

has Mellin transform: 

r(a + s)r(a + 6) „ 
r(a)r(a + o + s) 

which allows to transform 1)541) into 

E (^) = ( V) " E ^/4,l/4) S E ^5/4) S , K S > -3/4 , 

and then A l = ^ $3/4,1/4 ■ 75/4 or , coming back to M BST 



M BST (m law ^ (Aj/4>l/4) -V2 . (75/4) -l/2 (5g) 

where the two variables in the right-hand side are independent. 
Remark 6.1 

The limit martingale connection (3) rewritten with (|5U|) and (|55|) 
gives 

(73/2) 2 '= 4/3 3/4 ,l/4-7l-75/4- (56) 

This identity can be viewed as a consequence of the two following 
identities in law that can be found in Chaumont-Yor [_]. The first 
one is, for any parameter a 

( 7a ) 2 = 47a .70+1, (57) 

which is the stochastic version of (|53|) . and the second one (coming 
from the classical beta-gamma algebra) is 

7f '= %i-f -7i , (58) 
holding for a < 2. Both give 

(7a) 2 = 40» 1 _a.71.7ai! . (59) 

v / 2 2 2 

Identity (|56|) is obtained taking a = 3/2. 
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